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Abstract 

This paper presents a detailed study on analytical solutions to a general nonlin¬ 
ear boundary-value problem in finite deformation theory. Based on canonical duality 
theory and the associated pure complementary energy principle in nonlinear elasticity 
proposed by Gao in 1999, we show that the general nonlinear partial differential equa¬ 
tion for deformation is actually equivalent to an algebraic (tensor) equation in stress 
space. For St Venant-Kirchhoff materials, this coupled cubic algebraic equation can be 
solved principally to obtain all possible solutions. Our results show that for any given 
external source field such that the statically admissible first Piola-Kirchhoff stress field 
has no-zero eigenvalues, the problem has a unique global minimal solution, which is 
corresponding to a positive-definite second Piola-Kirchhoff stress S, and at most eight 
local solutions corresponding to negative-definite S. Additionally, the problem could 
have 15 unstable solutions corresponding to indefinite S. This paper demonstrates 
that the canonical duality theory and the pure complementary energy principle play 
fundamental roles in nonconvex analysis and finite deformation theory. 

Keywords: Nonlinear elasticity, finite deformation theory, nonlinear partial differential 
equation, analytical solutions, canonical duality theory 

1 Nonconvex Variational Problem and Motivation 

A large class of finite deformation problems in nonlinear elasticity can be formulated on the 
basis of a variational principle (V) in which it is required to minimize certain nonconvex 
potential energy. Typically, this takes the form 

{V) : min |n(x) = [ lT(Vx)dO-b [ (j){x)pdn- [ x-tdr|, (1) 

L Jn JQ Jrt J 
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where x represents the deformation field (a bijection), W(F) is the strain energy per unit 
reference volume, which is a nonlinear differentiable function of the deformation gradient 
F = Vx, and V is the gradient operator in a simply-connected domain (the reference 
configuration of the body) C with boundary 917 = F = Ft U F^ such that F^ n F^ = 0. 
Each material point in 17 is labeled by its position vector X and the corresponding point 
in the deformed configuration is denoted by x (= x(X)). The body force f (per unit mass) 
is taken to be conservative with potential so that f = —gradi;A, and p is the reference 
mass density. On the part F^ of the boundary the surface traction t is prescribed to be 
of dead-load type, while on F^ the deformation x is given. The notation Xa identifies a 
kinematically admissible space of deformations %, defined by 

^a = {xG W^’P(f7;M3) I Wx^J^a, X = Xq on F^}, (2) 

where is the Sobolev space, i.e. a function space in which both x a.nd its weak 

derivative Vx have a finite LP{VL) norm. = {F G T^(17;R^^^)| detF > 0} denotes the 
admissible deformation gradient space with p > 1. Clearly, solutions % G Ta of the problem 
{V) are not necessarily to be smooth. 

The criticality condition (5Fl(x) = 0 leads to a mixed boundary-value problem (BVP), 
namely 

f V • [VFkF(Vx)] + pf = 0 in 

(BVP) : (3) 

[ n • [VfPF(Vx)] = t on Ft, 

where VfVF(Vx) = dW{F)/dY (in components dW/dPia)-, n is the unit outward normal 
to Ft and, in component form, we adopt the conventions V • r = {dTia/dX^} and r • n = 
{Tiana}- Note that V-r is defined in the weak sense where Vx is discontinuous. In general, 
it is rarely possible to solve this nonlinear boundary-value problem by use of direct methods. 
Indeed, the strain energy W{F) is a nonconvex function of F, the problems {V) and (BVP) 
are not equivalent, and {BVP) may possess multiple solutions. Identification of the global 
minimizer of the variational problem {V) is a fundamentally difficult task in nonconvex 
analysis. From the point of view of numerical analysis, any numerical discretization of the 
problem {V) leads to a nonconvex minimization problem, and it is well known in global 
optimization theory that most nonconvex minimization problems are NP-hard [14, 15, 16]. 

Duality principles play fundamental roles in sciences and engineering, especially in con¬ 
tinuum mechanics and variational analysis. For linear elasticity, since the stored strain 
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energy W is a convex function of the (infinitesimal) strain tensor, it is well-known that each 
potential variational (primal) problem is linked a unique equivalent (dual) complementary 
variational problem via the conventional Legendre transformation. This one-to-one du¬ 
ality relation is also known as the complementary variational principle, which has been 
well-studied with extensive applications in both mathematical physics and engineering me¬ 
chanics (see Arthurs, Nobel-Sewell, Oden-Reddy, Tabarrok-Rimrott, etc). 

In finite deformation theory, if the stored energy density IT(F) is a strictly convex 
function of the deformation gradient tensor F over the field fl, then the first Piola-Kirchhoff 
stress tensor can be uniquely determined by r = VIT(F) and the complementary energy 
density W* can be obtained explicitly by the Legendre transformation: 

IT*(t) = {F:t-IT(F)| t = VIT(F)}, (4) 

where F : r is defined as tr(F • r"*") and signifies the transpose. In this case, the comple¬ 
mentary variational problem can be defined as 

min •|n‘^(T) = J W*{T)dQ. — J Xq ' ' J^dF , (5) 

where Ta is the statically admissible space defined by 

Ta = {t ^ C^{Vt) I V • T -|- pf = 0 in Li, T • n = t on Ft} , (6) 

where q is the conjugate number ofp, i.e. it is given by 1/p+^/q = 1. This complementary 
variational problem was first studied by Levinson [31]. The well-known Levinson principle 
states that if r is a solution of the complementary variational problem (5), then the defor¬ 
mation field X defined through the inverse constitutive law F(x) = VIT*(t) is a solution 
of the potential variational problem (1) and the complementarity condition 

n(x) + n'=(T) = 0 

holds. This principle can be proved easily by using the traditional Lagrangian duality theory 
(see Gao, 2000). 

The Levinson principle is simply the counterpart in finite deformation theory of the 
complementary variational principle in linear elasticity. In finite deformation theory, the 
stored strain energy IT(F) is in general nonconvex such that the stress-deformation relation 
T = VIT(F) is not uniquely invertible [36, 37] and the complementary energy function W* 
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cannot be defined explicitly via the Legendre transformation. Although by the Fenchel 
transformation 

VL“(t) = max{F:T - bF(F)}, 

F 

the Fenchel-Moreau type dual problem can be formulated in the form of 

m^ |n**(r) = J Xo '~ J fF**(T)dn|', (7) 

the nonconvexity of W leads only to the so-called weak duality theorem 

min nlv) > maxlf^fr) 

XeAfa reTa 

due to the Fenchel-Young inequality VF(F) > F : r — VF^(t). In nonconvex analysis, the 
nonzero 9 = n(x) — max-reT; n^(T) > 0 is called the duality gap. This duality gap 

shows that the well-developed Fenchel-Moreau duality theory can be used to solve mainly 
convex problems. 

In finite deformation theory, the well-known Hellinger-Reissner principle [25, 38] and 
the Fraeijs de Veubeke principle [46] hold for both convex and nonconvex problems. How¬ 
ever, these principles are not considered as pure complementary variational prineiples since 
the Hellinger-Reissner principle involves both the displacement field and the second Piola- 
Kirchhoff stress tensor; while the Fraeijs de Veubeke principle has both the rotation tensor 
and the first Piola-Kirchhoff stress as its variational arguments. The existence of a pure 
complementary variational principle in general finite deformation theory has been discussed 
by many researchers over several decades (see, for example, [28, 29, 30, 32, 35, 36, 37]). 
Moreover, since the extremality condition in nonconvex variational analysis and global op¬ 
timization is fundamentally difficult to resolve, none of the classical complementary-dual 
variational principles in finite deformation theory can be used for reliable numerical com¬ 
putations. 

Canonical duality theory provides a potentially useful methodology for solving a large 
class of nonconvex problems in complex systems. This theory consists mainly of (1) a 
eanonical dual transformation, which can be used to formulate perfect dual problems in 
nonconvex systems; (2) a complementary-dual variational prineiple, which allows a unified 
analytical solution form in terms of the canonical dual solutions; (3) a triality theory, which 
provides sufficient criteria for identifying both global and local extrema. The original idea 
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of the canonical dual transformation was introduced by Gao and Strang [23] in finite defor¬ 
mation systems. In order to recover the duality gap in nonconvex variational problems, they 
discovered a so-called complementary gap function, which leads to a complementary-dual 
variational principle in finite deformation mechanics. They proved that if this gap function 
is positive on a dual feasible space, the generalized Hellinger-Reissner energy is a saddle- 
functional. It turns out that this gap function provides a sufficient condition for global 
optimal solution in nonconvex variational problems. Seven years later, it was realized that 
the negative gap function can be used to identify local extrema. Therefore, a triality theory 
was first proposed in post-buckling problems of a large deformation beam model [7], and a 
pure complementary energy principle was eventually obtained in [9]. This principle can be 
used to obtain a general analytical solution for 3-D large deformation elasto-plasticity [11]. 
It was shown by Gao and Ogden (see [19, 20]) that for one-dimensional nonlinear elastic¬ 
ity problems, both global and local minimal solutions are usually nonsmooth and can’t be 
obtained by any Newton type of numerical methods. For finite dimensional systems, the 
canonical duality theory has been successfully applied for solving a large class of challeng¬ 
ing problems in computational mechanics [4, 24, 43] and global optimization with extensive 
applications in computational biology [47], chaotic dynamical systems [33, 41], discrete and 
network optimization [18, 21, 40, 42]. 

The purpose of this paper is to illustrate the application of the pure complementary 
variational principle in combination with triality theory by solving a general nonconvex 
variational problem governed by St Venant-Kirchhoff material. The paper is organized 
as follows. Section 2 presents a brief review on the canonical duality theory in nonlinear 
elasticity. Some fundamental issues in nonlinear elasticity are addressed, including the 
reasons why the Legendre-Hadamard condition provides only necessary condition for local 
minima, how the Gao-Strang gap function and the triality theory can be used to identify 
both global and local extremal solutions. In Section 3 we show that for the St Venant- 
Kirchhoff materials, the pure complementary variational problem can be solved principally 
to obtain all possible solutions. Some concluding remarks are contained in Section 4. 
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2 Canonical Duality Theory and Complementary Variational 
Principle 

It is known that the stored-energy function Vh ; —)■ M must obey certain physical laws and 
requirements in continuum mechanics, such as the principle of material frame-indifference 
[45], which lay a mathematical foundation for the canonical duality theory. Let 

SO(3) = {Q G qT ^ Q~i Q = 1} (8) 

be the special orthogonal group. 

Definition 2.1 (Objectivity and Isotropy [12]) 

(Dl) Objective Set and Objective Function: A subset Ta C is said to be 

objective if for every F G Ta and every Q G SO(3), QF G Fa- A scalar¬ 
valued function W : ^ IK is said to be objective if its domain is objective 

and 

IF(QF) = W{F) VF G Fa, VQ G SO(3). (9) 

(D2) Isotropic Set and Isotropic Function: A subset Fa C is said to be 

isotropic if for every F G Fa and every Q G SO(3), FQ G Fa- A scalar¬ 

valued function W ; Ja —^ IK is said to be isotropic if its domain is isotropic 
and 

IF(FQ) = W{F) VF G Fa, VQ G SO(3). (10) 

The objectivity implies that the constitutive law of material is independent with the 
observer (coordinate free). While the isotropy means that the material possesses certain 
symmetry. Generally speaking, the deformation gradient F is a two-point tensor, which 
is not considered as a strain measure. The right Cauchy-Green tensor C = F^F is a 
(Lagrange type) strain measure which is objective (rotation free), i.e., 

C(QF) = (QF)'^(QF) = F^Q^QF = C(F) VQ G SO(3). 

Dually, the left Cauchy-Green tensor B = F^F is an isotropic function of F. 

In continuum mechanics, the objectivity is also known as the principle of frame-indifference- 
According to P.G. Giarlet, the stored energy function of a hyper-elastic material is objective 
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if and only if there exists a function U{C) such that WiY) = [/(C(F)) (see Theorem 4.2-1 
in [5]). This principle lays a foundation for the canonical duality theory. 

Indeed, the canonical dual transformation was developed from the concept of the objec¬ 
tivity. The key step of this transformation is the introduction of a geometrically admissible 
strain measure ^ = A(x) ; Tjj —)■ C and the canonical function : Tq ^ M such 
that the nonconvex stored energy 1T(F) can be written in the canonical form IT(Vx) = 
U{A(x))- According to [12], a convex differentiable real-valued function U{^) is said to 
be canonical on its domain £a if the duality relation = Vt/(^) : £a ^ £a is invertible 
such that the conjugate function U*{^*) of U{^) can be defined uniquely by the Legendre 
transformation 

u*{e) = {^: r - um r = vt/(^) g t4. (h) 

By the theory of convex analysis, it is easy to prove that the following canonical duality 
relations hold on £a x £* 

^* = VU{^) ^ ^ = VU*{C) ^ + (12) 

and the pair (^,^*) is called the canonical dual pair on £aX £*. 

Thus, on replacing VL(Vx) in the total potential energy n(x) by its canonical form 
= U{A{x)), and we take the body force to be a constant, so that (pix) = the 

minimal potential energy variational problem (1) can be written in the following canonical 
form 

(V) : min |n(x) = /'[17(A(x))-/9X-f]dn- / x-tdr|. (13) 

Furthermore, in terms of and by the Fenchel-Young equality 

U(A(x)) = A(x):^-U*(^), 

the so-called total complementary energy funetional [23] H : Ajj x T* —)■ M can be written, 
in the present context, as 

S(X,<?) =/ [A(x) :<?-!/*(<?)- / X • tdF. (14) 

Jn JFt 

For a given statically admissible field r G Ta, this total complementary functional can be 
written in the following form 

St(x,<?)=/ Xo • -ndF-F / [A(x);<J-!/*(<?)-(Vx) ; t] dfl. (15) 

Jr^ Jfi 
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( 16 ) 


For a given G the canonical dual functional is then defined by 

I = 0 } = - / U*{^)dn, 

Jn 

where is dehned by the so-called A-conjugate transformation [12, 16] 

=sta|^[A(x):<j-px-f]dri-^ X-tdr I % G , (17) 

with sta indicating the stationary value at fixed <;££*. In terms of r £ Ta, we have the 
following form 

= Xo • • ndr-Fsta|^[A(x):<j - (Vx) : rjdfl I X G A^j . (18) 

In finite deformation theory, 

n^(<j) = F^(<j) - [ u*{^)dn (19) 

Jn 

is also called the pure complementary energy functional, which was first proposed in [9]. 

Theorem 2.1 (Complementary-Dual Variational Principle [11]) For a given stati- 
eally admissible field r G Ta, the following statements are equivalent: 

F (X)^) is a eritieal point o/H' 7 -(X)<?)/ 

2. X is a eritieal point ofIl{x); 

3. ^ is a eritieal point o/n^(<j). 

Moreover, we have 

n(x) = s(A, ?) = Sr (x, ^) = n^(?)- (20) 

This theorem shows that to find a critical solution to the nonconvex total potential 
n(x) is equivalent to find a critical point of its canonical dual function nr(<?). For a 
given T G Ta, different choice of the geometrical measure A(x) will leads to different, but 
equivalent, nr(<?) on a subset 5a C f*. 

In finite deformation theory, the canonical duality relation is also known as the Hill 
work conjugate and the canonical function U{$,) is called strain energy-density. According 
to Hill, for a given hyper-elastic material, there exist a class of strain measures ^ and the 



associated canonical functions U (^) such that the associated stress can by defined uniquely 
by the canonical duality relation = V?7(^). There are many canonical strain measures 
in finite elasticity and many of these strain measures belong to the well-known Hill-Seth 
strain family 

eW = 1[c'^-i], 

Irj 

where I is an identity tensor in and 77 is a real number. 

Canonical duality theory and pure complementary energy principle for general strain 
measures have been studied in [12], In this paper, we consider only the Green-St Venant 
strain tensor simply denoted as E. In this case, the geometrical operator 

E = A(x) = ^[(Vx)^(Vx) -l]-.Xa^£a (21) 

is a quadratic operator and its domain can be defined by 

fa = {E e ^3X3^1 ^ (2E + I)^0}. (22) 

We assume that the associated strain energy density U (E) : fa —^ is convex such that 

the conjugate stress ^ of E, denoted by S, can be defined uniquely by the constitutive law 


S = V[/(E) : fa ^ fa*. 


(23) 


This associated stress S is the well-known second Piola-Kirchhoff stress, which is well- 
defined on f* = {S G fPAp-2)(^f7.]^3x3^| g _ s'r|_ jji case, the pure complementary 
energy 11^ has the form of 

'1 


tr(T-S-^-T + S) + [/*(S) 


dn. 


(24) 


(25) 


n^(S)= [ Xo-T-ndT- [ 

which is well-dehned on the canonical dual space 

= {S G f*| tr(T • S-^ • r) G Vr G Ta}. 

Therefore, the canonical dual problem is to hnd the critical point S G fa such that 

(P'^) : n^(S) = sta{n^(S)| S G fa}. (26) 

Theorem 2.2 (Analytical Solution Form [12]) For a given r G Ta, if S is a critical 
point o/n^(S), then along any path from Xq G to X. £ FI, the deformation defined by 

X= T-S-i-dX + Xo(Xo) (27) 

4Xo 
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is a critical solution to {V). Moreover, if 

Vx(t-S"^) = 0, (28) 

then X is a closed form solution to the boundary value problem (BVP) (3). 

The proof of this theorem can be found in [8, 9, 11]. In fact, the criticality condition 
(5n^(S) = 0 leads to the following dual tensor equation: 

S-[I + 2(V[/*(S))]-S = t^-t, (29) 

which is equivalent to 

V17*(S) = S-^)^T • - I) . 

This is actually the constitutive law E = A(x) = i[F^F — I] = V?7*(S) subjected to 
F = T • S“^. Therefore, if the compatibility condition V x F = 0, in index notation 

dFia _ dFip 
dXp ~ dX^ ’ 

holds, then F is the deformation gradient and % is a solution to {BVP). 

Remark 2.1 (PDF Algebraic Equation) Theorem 2.2 shows that by the pure com¬ 
plementary energy principle, the nonlinear partial differential equation {BVP) is equiva¬ 
lently converted to a canonical dual tensor equation (29), which can be solved to obtain the 
stress field S for certain materials. From the equation (29) we know that S = 0 i/r = 0. 
Therefore, although appears in n^(S), this pure complementary energy is well-defined 
on Sa. The equation (27) presents an analytical solution form to the boundary value prob¬ 
lem in terms of the canonical dual stress field S and the statically admissible t G Ta. Of 
course, this is purely formal and in general it is not easy to obtain the solution for general 
practices unless the deformation compatibility condition (28) holds. 

It has been assumed here that the relation between S and E is invertible. This certainly 
holds in a neighborhood of the (stress-free) reference configuration since the canonical strain 
energy U (E) is convex in such a neighborhood. It is a reasonable assumption to extend 
this to a sufficiently large domain that includes deformations of practical interest. Finite 
element implementations of nonlinear elasticity are usually based on the variables S and 
E and the associated tangent tensor dS/dF = V^?7(E), which is assumed to be positive 
definite. It is always possible to select forms of the strain-energy function W such that this 
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is the case, although the possibility of its failure for particular materials is not in general 
ruled out. 

In terms of the deformation X ^ and the second Piola-Kirchhoff stress S £ £*, the 
total complementary functional H(x, S) can be written as 

Hr (X, S) = / [E(x) : S - C/*(S) - (Vx) : r] df^ + / Xo ' t • ndP (30) 

J £7 

which is actually the well-known Hellinger-Reissner energy if the first Piola-Kirchhoff stress 
is replaced by external force field. From the nonlinear canonical dual tensor equation (29) 
we know that for a given r G Ta, the pure complementary energy 11^(8) may have multiple 
critical points. In order to identify the global extremum, We need to introduce the following 
subspaces: 

5+ = {Sg 5,| S^O}, 5- = {SG5a| S^O}. (31) 

Theorem 2.3 Suppose for a given r G Ta, the pair (x^S) isolated critical point of 

Ht(X)S). //S G then X is a global minimizer o/n(x) on Xa if and only ifS is a global 
maximizer o/11^(8) on Sfi , i.e., 

n(x)=minn(x) ^ maxn^(S) = n^(S). (32) 

Xea'a se5+ 

// 8 G S~, then X is a local maximizer of n(x) if and only if S is a local maximizer of 
n^(8), i.e., on a neighborhood Xq x So C Xa x , 

n(X) = max n(x) maxn^(S) = n^(S). (33) 

S€<5o 

// 8 G S~ and VpIF(Vx) 0, then X is a local minimizer o/n(x). 

Remark 2.2 (The Complementary Gap Function and Triality Theory) 

Theorem 2.3 shows that the extremality of the primal solution x depends on its canonical 
dual solution 8. This result was first discovered by Gao and Strang in 1989 [23], i.e. they 
proved that X(S) is a global minimizer of n(x) if the complementary gap function satisfies 

Gapix, S) = / kvx)^(Vx) + I] : 8 dfl > 0 Vx G (34) 

Jn ^ 
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Since Gap(x, S) is quadratic in x> this gap function is positive for any given X ^ if 
S ^ 0. Replacing F = Vx & 2 / F = r • S~^, this gap function can be written as the so-called 
pure gap function 

GapixiS), S) = / l:tr{T ■ • T + S) d^l, (35) 

JQ ^ 

which is a main term in the pure complementary energy n^(S) in addition to U*{S). Com¬ 
paring n^(S) with nt*(T ) given by (7), we can understand that this gap function not only 
recovers the duality gap in the Fenchel-Moreau duality theory, but also provides a global 
extremality condition for nonconvex variational problem iV). 

To see this in detail, let us consider the canonical transformation W{¥) = C/(E(F)). By 
chain rule we have 


d‘^W{Y) 

dFfdFf^ 


3 

^ F^gHeapuFi, 
e,u=i 


(36) 


where H = {HoafSu} = V^C/(E). By the convexity of the canonical function UfE), we have 
H 0. Therefore, if S = {5'q,^} G the Flessian V^VF(F) 0 and, by Gao and Strang 
[23], the associated deformation field X is a global minimizer ofIl{x)- The statement (32) 
shows that the nonconvex minimization problem i[P) is equivalent to a concave maximization 
dual problem over a convex space Sfi, i.e.. 


max{n^(S)| S€5+}, 


(37) 


which is much easier than the nonconvex primal problem {V). The global optimality condi¬ 
tion S G 5^ is a strong case of Gao and Strang’s positive gap function (34)- 

Subsequently, in a study of post-buckling analysis for a nonlinear beam theory, it was 
found that if the dual solution S is negative definite in the domain ^l, the solution X could 
be either a local minimizer or a local maximizer of the total potential energy. To see this, 
we substitutive F = r • into (36) to obtain 

= E ^!lS3Wc0.S-iri (38) 

which shows that even if S ^ 0, the Flessian matrix V^VF(F) could be either positive or 
negative definite, depending on the eigenvalues of S G Sf. Thus, in addition to the double¬ 
max duality (33), we have the so-called double-min duality 

n(x) = min n(x) ^ min n^(s) = n^(s), (39) 
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which holds under certain condition (see [15]). For this reason, a so-called triality theory 
was proposed first in post-buckling analysis of a large deformed beam model [7], and then 
in general nonconvex mechanics [11, 12]. This triality theory reveals an important fact 
in nonconvex analysis, i.e. for a given statically admissible field r G Ta, if the canonical 
dual equation (29) has multiple solutions {S^} in a subset Qq C ri, then the boundary 
value problem (BVP) could have an infinite number of solutions {x;j(X)} in Q. The well- 
known Legendre-Hadamard (L-H) condition is only a necessary condition for a local minimal 
solution, while the triality theory can identify not only the global minimizers, but also both 
local minimizers and local maximizers. It is known that an elliptic equation is corresponding 
to a convex variational problem. If the boundary-value problem (3) has multiple solutions 
{X;j(X)} at one material point X G the total potential n(%) is not convex and the 
operator A{x) = V • [VfM^(Vx)] may not be elliptic at X. G H even if the L-H condition 
holds at certain ^^(X). 

The pure complementary energy principle and triality theory play a fundamental role 
not only in nonconvex analysis, but also in computational science and global optimization 
(see [15, 17, 18, 22]). 


3 Application to St Venant-Kirchhoff Material 


For St. Venant-Kirchhoff material, the canonical energy function t/(E) has the most simple 
form: 

[/(E) = ;utr(E2) + -A(trE)2. (40) 

2 

The second Piola-Kirchhoff stress depends linearly on the Green-St Venant strain via the 
Hooke’s law: 

S = V[/(E) = 2/xE + A(trE)I = H : E, (41) 


where H is the Hooke tensor for St Venant-Kirchhoff material. The complementary energy 
is 


[/*(S) 


J_tr('s2i_ 

Ap. ^ ’ 4//(3A + 2/i) 


(trS)' 


(42) 


and hence 

E = Vf/-(S) = Ts-^-i^(t..S)I.H--S. (43) 
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By the canonical dual tensor equation (29), we have 


+ 2S(H-^ : S)S = + -S^ - / (trS)S^ = r'^r. (44) 

/i 2/ij 

The diagonalization of this tensor equation leads to the following coupled quebec nonlinear 
algebraic systems; 

For convenience, we make the following substitutions in (45): 

Si = rf = i = 1, 2,3, 

and k = <1/3 (due to jj, > 0). So, the system (45) can be written as follows 

/c(?i + ?2 + ?3)?f = CTi, i = 1,2, 3. (46) 

3.1 Auxiliary Equation 

In this section we will study solutions of the following equation: 

G(?,g,fT) = ?^ + (1 -/cg)?^ - o-= 0, (47) 

where (T>0, 0</c< and q is an arbitrary real number. Also, since a > 0, we can 
assume that ? / 0. 

Since the parameter q in this section is assumed to be independent on the following 
results are similar to one-dimensional nonlinear elasticity problems studied by Gao [12, 13], 
Gao and Ogden [19]. 

Lemma 3.1 If o,re solutions of the equations G{(^,q,ai) = 0, G(?,g,CJ 2 ) = 0, 

G(?, g, as) = 0 correspondingly, and + ?2 + ?3 = Q, then ?i, ? 2 , ?3 satisfy (46)- 

Proof. Obvious. □ 


Lemma 3.2 Equation (47) has exactly one positive solution. It has negative solutions iff 





There is only one negative solution if and only if q 


1 ( 1-3 


\/i)* 
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Proof. To check that there is exactly one positive root one can apply the Descartes’ rule of 
signs. To prove the rest, let’s fix q, a and notice that G(?, g, a) = 0 has negative solutions iff 
it has at least two different solutions. This will happen iff the values of the function at local 
minimum and maximum have different signs. The extremums of G are at ?o = —1(1 — kq) 
and 0. Since the value of G at 0 is —a < 0, we find when G{<;o,q,a) > 0. Solving this 
inequality we get <? < -^(1 — 3^/^). 

Corollary 3.1 The equation G(?,0, ct) = 0 has negative solution(s) iff 

4 

a < —. 

“ 27 

Proof. Apply Lemma 3.2 to g = 0. □ 

Lemma 3.3 Let’s fix a > 0 and assume that go,qo satisfy (47), o.nd <Jo / 0, ?o / 

Then there exists a unique continuously differentiable function g{q), such that ?(go) = ^o, 
<;{q) and q both satisfy (47) and 

dg kifi 
dq <^3 + 2a 

Moreover, there are three possibilities (“branches”) forg{q): 

(a) If go G (—oo,—v^^), then the range of g{q) is the domain is 

(—oo, -^(1 — and g{q) is monotonically increasing. 

(b) If go £ then the range of g{q) is the domain is 

(—oo, -^(1 — ^f/j))) , and g{q) is monotonically decreasing. 

(c) If go G (0,+oo), then the range of g{q) is (0,+oo), the domain is (—oo,+oo), 
and g{q) is monotonically increasing. 


Proof. Let’s fix a and find q from (47) 




0-^ + 0 - cr 


kg"^ 


Since, ^ and a > 0 it is obvious that q{g) is monotonically increasing in the 

intervals o G (—oo, — v^) and o G (0, +oo) and is monotonically decreasing in the interval 
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G (—0). The corresponding intervals for q are (— 00 , ^(1 — 3^/j)), (— 00 ,+ 00 ), and 
(— 00 , ^(1 — Also, one can easily check that 

d? _ 

dq~ + 2a' 

Thus, the lemma is proved. □ 

Definition 3.1 The three branches of <;{q,a) (a is fixed) described in Lemma 3.3 will be 
denoted as follows: 

(a) <^^{q,a) is a positive branch with the domain (— 00 ,+ 00 ) and range (0,+ 00 ); 

(b) cr) < cr) are two negative branches with the domain (— 00 , -^(1 — 3^/|")) and 
ranges (—oo,—\'^) and (—0) correspondingly. (Note that Corollary 3.1 implies 
that i(l -3^) < 0 ; 

Definition 3.2 Let’s introduce the following notations: 

(a) e(g, a) = fi{q, a) - f = 1, 2,3; 

(b) F’-’^’’^{q,ai,a2,a3) = ^{q,ai) + {q,a2) +g'^{q,a3), i,j,k = 1,2,3. 


Lemma 3.4 The following statements are true: 
(a) For i = 1,2, 3 


= - 


(1 - Skfifiq, afi + fi{q, afi - a 


and 


3kC{q, afi 
dC {l-3k)fi{q,afi + 2a 


dq 3(<?*(g, cj)3 + 2cr) 

(b) ?^(0,cr) = ?^(0,fj) > 0, ?^(0,cj) = ?^(0 ,(t) < 0, and ?^(0,a) = g^{0,a) < 0. 

(c) For a fixed a, g^{q,a) is monotonically decreasing in q and 

limg^+00 = -00 • 


(d) For a fixed a, 


lim <i^{q,a) = +00 and lim ^'^{qfi) = + 00 . 

q—^—oo q—^—00 
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(e) For fixed ai,a 2 ,o's, each o/F*’-^’^(g, cji, < 72 , 1 T 3 ), i,j,k = 1,2,3, is eontinuous. Moreover, 
cji,( 72 ,< 73 ) is monotonieally decreasing in q. 

Proof. To check (a), first substitute q{g) = into ?'{q,cr) = ?*(?, 7 ) — |, i = 1,2,3. 

Expression for ^ can be obtained either by direct differentiation of the previously obtained 
expression for 7 ) or subtracting ^ from ^ 

(b) is obvious. 

To prove (c), recall, that k < ^, and use formulas from (a). To prove (d), recall, that k < ^, 
and use the first formula from (a), (e) immediately follows from (a) and (c). □ 

Lemma 3.5 Solutions, ?^(0, 7 ), ?^(0, 7 ), ?^(0, 7 ), of the equation G{q, 0, 7 ) = — 7 = 

0, 0 < 7 < ^, enjoy the following properties: 

(a) If a = -^ the solutions are ?^(0, ?^(0, = ?^(0, | 

(b) IfO < 7i < 72 < 0, then 

0 < ?H0,7i) < ?^(0,72) < ^ 

and 

-1 < ?^(0,7l) < ?^(0,72) < < ?^(0>O-2) < ?^(0,7l) < 0 

(c) ?H0,7) + ?^(0,7) < 0 
Proof. (a) can be checked directly. 

To prove (b), one can either apply the implicit function theorem to H{<^, 7 ) = G{g, 0, 7 ) = 0. 
Or, less formally, draw the graph of y = and observe what happens to its roots 

when the graph is shifted upward until it becomes y = 

(c) Obviously, <^^(0, 7 ) + ?^(0, 7 ) + <^^(0, 7 ) = -1. So, ?^(0, 7 ) + <^^(0, 7 ) = -l-<^^(0,7) < 0, 
since ?^(0, 7 ) > —1. 

3.2 Solutions of the St. Venant-Kirchhoff Material 

We are now ready to present our main results. 

Theorem 3.1 For any given foree field f : fl —)■ and the surface traction t : T* —)■ 

such that the statieally admissible stress t G Ta has no zero eigenvalues almost ever where 
in Q, the canonieal dual problem {V^) has a unique positive critieal solution S G 
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Proof. We need to prove that for arbitrarily given cJi, 1 T 2 , (T 3 > 0, the system of equations 
(46) has a unique positive solution (?i,? 2 ,? 3 ), such that all q > 0, i = 1,2,3. From Lemma 
3.4(b), it follows that 

F^’^’^( 0 ,O-i,cr 2 ,cr 3 ) > 0 . 

From Lemma 3.4(c), it follows that for some qi > 0, large enough, 

F^’^’^(gi,cri,o-2,o-3) < 0 . 

Therefore, since T^44 jg continuous and monotonically decreasing in q (Lemma 3.4(e)), 
there exists a unique qo, Q < qo < qi, such that 

i.e. 

?H90,O-i) +?H90,0-2) +?^(^0,O-3) = ^ 0 - 

So, from Lemma 3.1 it follows that ?^((7o,<7i), <;^{qo,a 2 ), ‘i^{qo,o' 3 ) form a positive solution 
of (46), which are eigenvalues of the second Piola-Kirchhoff stress S. Therefore, Problem 
(V^) has a unique global maximizer S G 5^. □ 

Theorem 3.2 For any given force field f : 11 —)■ and the surface traction t : F^ —)■ 

such that the eigenvalues of the statically admissible stress tensor function r G Ta satisfy 
0 < (Ti,iT2,(T3 < the total complementary energy n^(S) has eight negative solutions 
Sfc G 5 -, A: = 1,... ,8. 

Proof. We need to prove that for arbitrarily given 0 < cJi, < 72 , CT 3 < ^, the system of 
equations (46) has 8 solutions (<?i,? 2 ,? 3 ), such that all q < 0, i = 1,2,3. From Corollary 
3.1 it follows that each of the equations G(<j, 0, CTj), has two negative solutions, ?^(0, cTj) > 
<?^(0, cjj), i = 1, 2, 3. From Lemma 3.4(b), it follows that for i,j, k = 2,3 

F*’^’'=(0,ai,c72,U3) < 0 . 

From Lemma 3.4(d) it follows that there exists < 0 such that 

F"’^’^{qi,ai,a2,a3) > 0 . 

Therefore, since F'^dF jg continuous in q (Lemma 3.4(e)), there exists qo, 0 > qo > qi, such 
that 

F"’^’^{qo,ai,a2,(T3) = 0 , 
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i.e. 


<?*( 9 o,o-i) +?^( 9 o,o- 2 ) + ^^{qo,(^3) = qo- 

So, from Lemma 3.1 it follows that ?*(go,ci), (fZO)'^ 2 ), ^^{qo,<^ 3 ) form a negative solution 
of (46). Since, each of i,j, k can be chosen independently from the set { 2 , 3}, we have total 
8 different negative solutions. □ 

Theorem 3.3 For any given force field f : —)■ and the surface traction t : —)■ 

such that the eigenvalues of the statically admissible stress tensor function r G Ta satisfy 
0 < cti,cj2,cj3 < the total complementary energy n^(S) has at least 15 mixed stationary 
points, i.e., some eigenvalues i = 1,2,3, 0 /S are positive, some are negative. 

Proof. Each of the equations G(?, 0,cjj), has one positive and two negative solutions: 
Applying Lemma 3.5 it is easy to check that 

(1) for i,j = 2,3, 

F^’*’^(0,cri,(T2,(T3) < 0,F*’^’^(0,CJi,cr2,CJ3) < 0 

( 2 ) < 0 , ^^’^’^( 0 ,ui,0-2,0-3) < 0 , ui,0-2,0-3) < 0 . 

(3) Fhl,2(0,fji,O-2,O-3) < 0, F^4,3(0,fji,O-2,O-3) < 0, F^’^4(0,fTi,O-2,O-3) < 0, 

F^’^’^(0,o-i,o-2,o-3) < 0 . 

For each of these 15 combinations, there exists < 0 such that cji, < 72 , < 73 ) > 

0. 

Therefore, since is continuous in q (Lemma 3.4(e)), there exists q^, 0 > go > such 
that 

F“’'’'^(go,/i,/2,/3) = 0 

that is 

?“(^o, o-i) + ?^(^o, 0 - 2 ) + ?‘'(^o, 0 - 3 ) = qo 

So, from Lemma 3.1 it follows that <j“(go,<7i), <^^{qo,o' 2 ), ‘i^{qo,o' 3 ) form a mixed solution of 
(46). 

Obviously, these 15 combinations result in different mixed stationary points of 11^ (S). □ 
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4 Conclusions 


We have illustrated that by using the canonical duality theory, the nonconvex minimal 
potential problem (P) is canonically dual to a concave maximization problem in a convex 
stress space , which can be solved by well-developed numerical methods. By the pure 
complementary energy principle, the general nonlinear partial differential equation in non¬ 
linear elasticity is actually equivalent to an algebraic (tensor) equation, which can be solved 
for certain materials to obtain all possible stress solutions. Both global and local extremal 
solutions can be identified by the triality theory, while the Legendre-Hadamard condition is 
only necessary for local minimizers. Our results shows that for St. Venant-Kirchhoff mate¬ 
rial, the nonlinear boundary value problem could have 24 solutions at each material point, 
but only one global minimizer if the statically admissible stress r 7 ^ 0. It is important to 
have a detailed study on these solutions in the future. 
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